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The escape of an ensemble of particles from the Bunimovich stadium via a small hole has been 
studied numerically. The decay probability starts out exponentially but has an algebraic tail. The 
weight of the algebraic decay tends to zero for vanishing hole size. This behaviour is explained by the 
slow transport of the particles close to the marginally stable bouncing ball orbits. It is contrasted 
with the decay function of the corresponding quantum system. 
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I. INTRODUCTION 

The decay of a quantum mechanical system whose 
states may be treated statistically is — on the average 
over the states — not always exponential. It is alge- 
braic, if only a few decay channels are open. This has 
been demonstrated in Q). There, the quantum mechani- 
cal system was simulated experimentally by a quasi-two- 
dimensional microwave cavity. The classical analog to 
this is the motion of a particle in a billiard shaped like a 
stadium which was shown by Bunimovich |^| to be fully 
chaotic and especially ergodic. By the present note, we 
want to show that the algebraic decay observed in Q is a 
quantum mechanical feature that has no classical coun- 
terpart. To this end we have studied the escape probabil- 
ity P(t) via a small hole out of the classical Bunimovich 
stadium and find P(t) "almost exponential" . This means: 
Every P(t) has an algebraic tail for large t, but the decay 
functions approach exponential behaviour with decreas- 
ing size of the hole. 

The algebraic asymptotics of the decay function P(t) 
can be viewed as an example of the "third type of decay" 
mentioned in the introduction of H , where "particles ini- 
tialized in a chaotic region can stick a long time to the 
vicinity of the boundary of a regular region" . Although 
there is no region of regular motion in the phase space 
of the billiard there exist marginally stable orbits that 
cause this "sticking". A schematic model will let us - 
semiquantitatively — understand why the decay is nev- 
ertheless almost exponential. 

The results of the numerical simulation are presented 
in Sec. H. They are contrasted with the decay of the 
corresponding quantum system in Sec. III. The simple 
model illustrating the interplay between exponential and 
algebraic decay of the classical system is described in 
Sec. IV. Its predictions are compared with the present 
results in Sec. |v|. 



II. NUMERICAL EXPERIMENT 

As in [jjj, a quarter of the Bunimovich stadium has 
been considered in order to remove symmetry. Its shape 
is sketched in the insert of Fig. [|. All lengths are given in 
units of the radius of the circular part of the boundary. 
The shape parameter 7 (see Fig. l]) was 7 = 1.8. An 
escape hole of size A has been assumed — as indicated 
in Fig. [l] — in the upper half of the small straight piece 
of the boundary. For an ensemble of 10 6 particles, initial 
conditions were chosen at random. The distribution was 
a constant times dxdydcj), where x and y are Cartesian 
coordinates of the position inside the stadium and <fi is an 
angle characterising the direction of motion. The orbit of 
every particle was followed numerically until it escaped 
via the hole out of the stadium. The orbits were, however, 
not followed beyond 10 5 collisions between the particle 
and the boundary. For a hole of size A = 0.05 the result 



of this numerical experiment is given in Fig. @ by the 
histogram of the probability density P{L) for the particle 
to escape after an orbit of length L. This will be called 
data in the sequel. Since the velocity of the particle has 
constant modulus, we identify L with time. 

One finds P(L) exponential for sufficiently small L. 
Later P(L) turns into a more slowly decaying function. 
This behaviour is typical and it occurs wherever the hole 
is positioned. The probability density P(L) can be rep- 
resented by the combination of exponential and algebraic 
functions 



P(L) = EXexp (-XL) + Aa(fi - 1)(1 + aL)~ 



(1) 



where E + A = 1 for a normalized P. Then A measures 
the weight of the algebraic term. The full curve in Fig. |l| 
is a fit to the data. The parameters a, /3, A and A were 
searched for. This procedure was followed in 15 "exper- 
iments" with holes in the range of 0.25 > A > 0.0025. 
The results are partly reproduced in Tab. || together 
with the values of the normalized x 2 - Although only 
the 7 = 1.8 stadium is discussed and analyzed below, we 
display in Tab. [jjalso some results pertaining to the 7 = 1 
stadium in order to show that it behaves quite similarly. 

Table | and Fig. || show that the weight A of the alge- 
braic decay approaches zero for A — > 0, i.e. the decay of 
the stadium becomes exponential in the limit of a van- 
ishing size of the hole. We call this behaviour "almost 
exponential" decay. 

In the limit of A — » the decay constant approaches 
— see Tab. | — the value 
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given in Here, A c is the area of the billiard; the 
momentum of the particle has been set equal to unity. 
Equation (^) can be derived from ergodicity: Every point 
in phase space should be visited with equal probability 
independent of the elapsed time. From this one infers 
|| that the decay should be exponential with the decay 
constant Ao- 

We note that the dynamics of the stadium may also be 
described as a mapping which generates — from given 
coordinates of collision with the boundary — the coor- 
dinates of the next collision. If one uses the number of 
mappings applied to the initial distribution as time vari- 
able one again infers exponential decay, since the motion 
on the boundary is ergodic, too. The decay constant Vq 
is then simply the ratio of the size of the hole and the 
perimeter of the billiard. 
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(3) 



Here, we consider as boundary of the quarter stadium 
only the piece that it shares with the original full sta- 
dium. We disregard the "lines of cut" . Similarly a "colli- 
sion with the boundary" is defined as a collision with the 
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piece of the original boundary (one could as well include 
the full boundary of the quarter stadium into these def- 
initions but the present choice is more convenient). The 
ratio 

va_ = irjir + 4 7 ) 

An 2(tt + 2 7 ) 1 1 

then is the mean free path of the particle between two 
collisions with the boundary. We have verified numeri- 
cally that this relation is fulfilled (again in the limit of 
A — > 0). Equation fl) is a special case of a formula in 
the Appendix B of gf. 

The power j3 of the algebraic part of Eq. (0) is found 
to lie between 2 and 3. The fitted decay functions for 
7 = 1.8 with the parameters of Tab. | — which are use- 
ful representations of the data — are given in Fig. |^ as 
full lines. One notes that with A — > the onset of the 
algebraic decay is shifted to larger L and to decreasing 
probability level. 



it will decay essentially algebraically. Slight deviations 
from Gaussian distribution have been detected || and 
are due to the presence of bouncing ball orbits between 
the parallel straight sections of the stadium billiard. 

Hence, the essentially algebraic decay of the quantum 
system has no classical counterpart. The classical system 
with a small hole decays almost exponentially. 

What then is the origin of the algebraic tail in the 
classical decay functions P(L) ? — It is the fact that 
the motion in the stadium is ergodic only if it persists 
indefinitely If the motion terminates by the escape of 
the particle, it becomes apparent that ergodicity is not 
established in any finite time. The reason for this is the 
existence of marginally stable orbits that can almost trap 
the particle — the bouncing ball orbits. For details see 
the next section. 

It is then the same class of orbits — the marginally 
stable orbits — that add corrections to the essentially 
different behaviour of the classical and the quantum sys- 
tems. 



III. DECAY OF THE CORRESPONDING 
QUANTUM SYSTEM 

The almost exponential decay found above is in con- 
trast to the behaviour of the corresponding quantum sys- 
tem studied with the help of microwaves in (!]]. Suppose 
that the holes through which the system is coupled to the 
external world (the antennas) are small compared to the 
wave lengths occurring inside. Every hole may then be 
identified with one decay channel c. Let (r c ) be the de- 
cay width into channel c averaged over the eigenstates of 
the stadium. We suppose that the decay widths into dif- 
ferent channels are uncorrelated. To avoid complications 
which are not instructive, we take (r c ) to be independent 
of c. The probability P q that the quantum system decays 
at time t after its formation is then 



P q {t) ~ (1 + 2(T C )*)- 



(5) 



for t larger than the equilibration time of the system, see 
Sec. (6.3) of § or Eq. (6) of @]. Here, M is the num- 
ber of open channels (or holes). The algebraic form of 
Eq. (||) occurs because P q is an average over exponen- 
tially decaying resonances whose decay amplitudes have 
a Gaussian distribution. (See the introduction of 

The quantum system decays almost exponentially if M 
is large because Eq. (j^) can then be approximated by 



P q {t) ~exp [-(4 + M)(r c )t] 



(G) 



and the weight of the algebraic tail becomes negligible. 
For the channels to be statistically independent, any two 
antennas must, however, be separated by a typical wave 
length or more. 

We conclude that the quantum mechanical counterpart 
of a system with a small hole decays algebraically. If the 
statistics of the decay amplitudes is not exactly Gaussian 



IV. MODEL OF DELAYED CLASSICAL DECAY 

The argument justifying exponential decay with the 
decay constant Ao of Eq. (||) applies if ergodicity is es- 
tablished "sufficiently quickly". A close inspection |J 
reveals that the fluctuations of the frequency of the parti- 
cles' arrival at location of the hole must be small: Let At 
be the difference between two successive times of arrival 
at the hole. Then (At) 2 /At must be small compared to 
unity. Now, when there is a region of phase space into 
which the particles penetrate very slowly and — by the 
same token — in which they remain trapped for an ex- 
ceptionally long time once they are there, then the above 
fluctuations are large and the release of the trapped par- 
ticles will eventually dominate the decay process. Such 
regions of phase space have been described e.g. in jb] [l^] 
and the first Ref . of |2j : Close to the family of bouncing 
ball orbits and the "whispering gallery orbits" there are 
parts of phase space with volume > in which the parti- 
cles can be trapped for an arbitrarily long time. (A whis- 
pering gallery orbit is the motion of the particle along the 
boundary) . 

The arguments of |l(J show that the probability G(n) 
for the particle to be in an "almost bouncing ball or- 
bit" that will persist for more than n collisions with the 
boundary is 



G{n) 



(it + 2^)n ' 



n > 1. 



(7) 



An "almost bouncing ball orbit" allows for the angle be- 
tween orbit and boundary to be slightly different from 
7r/2. In this situation, n collisions amount to an orbit of 
length L — 2n. (See the definition of a "collision with 
the boundary" in Sec. |l|). Therefore 
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G(n(L)) = 



aoL' 



L > 1, 



(8) 



is the probability for the particle to be in an orbit that 
will persist over a length > L. Here, we have used 



ao 



TT + 2j 

2 7 2 ' 



(9) 



The probability to be in a whispering gallery orbit de- 
creases more strongly than L~ 1 and is therefore disre- 
garded in the further discussion. 

Hence, the probability that the particle is in an almost 
bouncing ball orbit that will persist for the length L is 
asymptotically 



9as(L) = ~G{n(L)) = 

aL aoL z 



(10) 



It was therefore anticipated in ||[14j that P(L) should 
asymptotically be ~ L~ 2 . There is some numerical evi- 
dence for this in [Q: Note that the function N(t) given 
there by closed circles in Fig. 2 is the present G(n(L)). 
The present data confirm this asymptotic behaviour of 
P(L) at least qualitatively. 

The following schematic model — inspired by the 
treatment of the Sinai billiard in jl3| — shows how the 
algebraically delayed decay comes about. Suppose that 
the phase space can be split into two parts C and C such 
that the decay happens in C, the delay in C. Once the 
particle is in C, it shall escape with probability u> or im- 
mediately go back to C with probability 1 — w. Consider 
a particle which may be anywhere at the time L = 0. 
Define g{L)dL as the probability distribution for its next 
transition from C into C to happen at time L. Then the 
probability density pi(L) to escape at time L after having 
made exactly one transition from C — > C is 



Pi(L) = ujg(L). 



(11) 



Obviously g as in Eq. ( ^0| ) is the asymptotic form of g. 
Let f(L)dL be the distribution of the time L between 
two successive transitions C — > C. Then the probabil- 
ity density P2{L) for escape after exactly two transitions 
C -► C is 



p 2 {L) = w(l - u)g <g> /, 



(12) 



where the operator ® denotes the convolution. For A; 
transitions C — > C one has 



Pk(L)=u(l-w)*- 1 g®(f®...f), 



(13) 



where g is folded with (k—2) convolutions of / with itself. 
The decay probability P(L) is the sum 



P(L) = £p fc (L). 



(14) 



k=l 



Taking the Laplace transform 



Pk(s) 



/•OO 

/ dLexp(-sL)p k (L) 
Jo 



(15) 



converts the convolution into a product of the Laplace 
transforms g and /. Therefore the transform P(s) of P 
is a geometric series in / which gives 



P(s)=ujg(s)[l-(1-Gj)f(s)} 



(16) 



The function g is by definition proportional to an in- 
tegral over /. We note that this is numerically exhibited 
by ijnj, where the function Pf pt corresponds to / and 
its asymptotic behaviour is found to be close to ~ L~ 3 . 
The relation between g and / is also just the relation 
between "transient chaos" and "chaotic scattering" es- 
tablished in jL2|. As the terms are used there, transient 
chaos corresponds to placing at time zero the particle at 
random in phase space; chaotic scattering corresponds to 
the situation in which the particle enters C at time zero. 
The results of |ll| are compatible with / ~ dg/dL as it 
should be. 

The normalization of / and g requires 



dL 



-fl(0)/(£) 



(17) 



and the definition of the Laplace transform then leads to 



/(*) =1-8 



5(0)' 



(18) 



Inversion of the Laplace transform (Eq. (|l6|)) then gives 
P(L) in the form of the Fourier integral 



P(L) = (2tt)- 1 ds'P(is')exp(is'L) 



(27T)" 



ds'- 



9(is') 



T i-l- ,~ r a exp(is'-L). (19) 
+ — —is'q (is') 



A similar formula was used in Eq. (5) of p5| to describe 
the exponential part only of the decay and to define the 
"relaxation time" within the closed stadium. It is shown 
here that Eq. ( Jig| ) accounts for the decay function alto- 
gether if the appropriate probability density g(L) is in- 
troduced. It must necessarily allow for large fluctuations 
of L. There exists of course a g(L) such that the observed 
decay function is exactly reproduced. The present model 
is schematic by a naive choice of g(L) leaving open the 
exact definition of the partition of phase space into C and 
C. We only know the asymptotic form of G, see Eq. (||). 
Let us — in the simplest possible way — assume that 



j(L) = ao(l + a L)~ 



(20) 



The decay function of eq. ( |19[ ) approaches an exponential 
when u> — > 0. The factor (1 — w)/wao in the denomina- 
tor of the integrand is then very large and this makes the 
functional form of g irrelevant. One can thus approxi- 
mately put g(is') w g(0) = 1 and the pole of P(is') at 
s' w iuja>o yields P{L) ~ exp (— ojo>oL). 
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The asymptotic behaviour of P(L) is given by setting 
the denominator of the integrand equal to unity. This is 
so because the omitted terms are of the type {is') n g{is') 
leading to derivatives of g(L) which decrease faster than 
g(L). Hence, P(L) — > g(L) for L — » oo. This means that 
the asymptotic behaviour is independent of the size of 
the escape hole! This result also implies P(L) ~ L~ 2 for 
large L as expected. Hence, the present model incorpo- 
rates the features of the data that we have termed almost 
exponential decay. 

By equating w«o with the decay constant Ao, Eqs. (|^) 
and (j^) determine the two parameters of the present 
model. This yields the theoretical decay functions given 
as dashed lines on Fig. ||| together with the data. 



V. DISCUSSION 

The numerical experiment yields an almost exponen- 
tial decay of the classical Bunimovich stadium. This 
means: The decay function P(L) starts out exponentially 
and later turns to an algebraic behaviour. This contrasts 
with the claim of |J] (see also the discussion in Jl(J). The 
weight of the algebraic decay, however, tends to zero with 
decreasing size of the escape hole. In this modified sense 
of exponential decay we agree with j| . 

Hillermeier et al. || show how the algebraic decay can 
be understood formally. The authors of [jl7 18 have 
pointed out that the slow transport of particles at the 
boundary of islands of regular motion can cause alge- 
braic decay. Yet there are no strictly stable orbits in the 
Bunimovich stadium [pj ; in this sense it is a system with 
"fully developed chaos" |Q. However, there is the family 
of marginally stable bouncing ball orbits which accord- 
ing to |L0| causes the algebraic damping of phase space 
correlations. It is responsible for the algebraic tail of the 
decay function P(L), too. More precisely, the bouncing 
ball orbits are a set of parabolic |2(|, non-isolated peri- 
odic orbits. A similar set of periodic orbits exists in the 
Sinai billiard and causes similar effects p3[ . 

The algebraic tail of P(L) and its suppression in the 
limit of vanishing size A of the hole — briefly the almost 
exponential decay — is scmiquantitatively reproduced by 
the model inspired by and described in Sec. IV. The 



]13| and described in Sec. 
model reproduces the exponential part of P(L). It should 
yield a lower limit to the algebraic part because there 
may be sources of delayed decay other than the bounc- 
ing ball orbits. The comparison between data and model 
in Fig. U shows this to be true. Furthermore the data 
and the results of the model should converge for L — > oo 
since the bouncing ball orbits cause the most pronounced 
delay. Convergence is indeed indicated by the curves for 
the largest holes although it is very slow; it takes many 
inverse decay constants Aq 1 before it is reached. One 
expects Q that P(L) tends — for L — > oo — towards 
a function g as ~ L~ 2 which is independent of the size of 
the hole and given by the properties of the bouncing ball 



orbits alone. The model of Sec. IV complies with this 
expectation. Again it is also compatible with the data: 
The experimental P(L) asymptotically behaves as L~ l3 
with [3 somewhat larger than 2 and thus indicate conver- 
gence towards g as (L). For the largest hole (A = 0.25) 
the convergence is essentially achieved within the range 
of L that we have studied. It is, however, achieved so 
slowly that for the other hole sizes A, even orbit lengths 
of L = 10 4 are insufficient to demonstrate it. The data 
are, however, fully consistent with convergence towards 

9as(L). 

There are numerous studies of algebraically delayed de- 
cay of correlations (as e.g. velocity correlations) in closed 
chaotic systems, see e.g. P, |l0| , ^l| , p2[ . The authors trace 
the delay back to the existence of marginally stable pe- 
riodic orbits. Anomalous di ffus ion is caused by the same 
type of orbits, see e.g. [E3 24 . Again these orbits are 



responsible for the algebraically delayed escape of parti- 
cles from fully chaotic systems. This emerges from the 
arguments of ]l4| and the numerical experiments of |l3| ] 
and the present note. Hence, all these phenomena are 
related to each other. By using the considerations of |lC| ] 
to define g as (L) we have directly linked the decay of cor- 
relations to the escape. 

Furthermore, we have pointed out that the algebraic 
decay of chaotic quantum systems is a consequence of 
wave mechanics and is not produced by the delayed decay 
of the classical counterpart. 
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TABLE I. Results of parametrizing the numerical experi- 
ments on the stadiums with 7 = 1.8 and 7 = 1. The quan- 
tities listed are those of Eqs. (^) and (Q) together with the 
normalized x 2 - va hies. 



y = 1.8 : 
A 


A 


a[HT 3 ] 


P 


A[10" 3 ] 


A [10- 3 ] 


x 2 


0.2500 
0.0500 
0.0100 
0.0025 


0.301 
0.181 
0.038 



34.00 
4.32 
4.71 

/ 


2.36 
2.89 
1.97 

/ 


34.93 
6.14 
1.22 
0.32 


30.77 
6.16 
1.23 
0.31 


0.99 
0.61 
0.50 
0.82 


7 = 1: 
A 


A 


a[KT 3 ] 


P 


A[10" 3 ] 


Ao[HT 3 ] 


x 2 


0.2500 
0.0500 
0.0100 
0.0025 


0.223 
0.128 
0.030 



48.66 
6.40 
2.09 

/ 


2.42 
2.96 
2.81 

/ 


49.93 
8.88 
1.79 
0.47 


44.54 
8.91 
1.78 
0.46 


1.38 
1.10 
1.03 
1.01 



FIG. 1. Decay probability of the stadium with 7 = 1.8 and 
an escape hole size of A = 0.05. The histogram is the result of 
the numerical experiment. The full line is its parametrization 
via Eq. (|l|). The parameters are given in Tab. (|j|). The dashed 
line is produced by the model of Sec. IV. 



FIG. 2. The weight A of the algebraic decay — see Eq. (|l]) 
— vs. the size A of the escape hole. The errors are statistical 
ones. They are given if larger than the size of the dots. 



FIG. 3. Decay probability P vs. orbit length (or time) 
L with hole size A as parameter. The full lines show the 
parametrized numerical experiment. The dashed lines result 
from the schematic model described in the main text. 
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